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Abstract 



This work is concerned with planar real analytic differential systems with 
an analytic inverse integrating factor defined in a neighborhood of a regular 
orbit. We show that the inverse integrating factor defines an ordinary dif- 
ferential equation for the transition map along the orbit. When the regular 
orbit is a limit cycle, we can determine its associated Poincare return map 
in terms of the inverse integrating factor. In particular, we show that the 
multiplicity of a limit cycle coincides with the vanishing multiplicity of an 
inverse integrating factor over it. We also apply this result to study the ho- 
moclinic loop bifurcation. We only consider homoclinic loops whose critical 
point is a hyperbolic saddle and whose Poincare return map is not the iden- 
tity. A local analysis of the inverse integrating factor in a neighborhood of 
the saddle allows us to determine the cyclicity of this polycycle in terms of 
the vanishing multiplicity of an inverse integrating factor over it. Our result 
also applies in the particular case in which the saddle of the homoclinic loop 
is linearizable, that is, the case in which a bound for the cyclicity of this 
graphic cannot be determined through an algebraic method. 

2000 AMS Subject Classification: 37G15, 37G20, 34C05 

Key words and phrases: inverse integrating factor, Poincare map, limit cycle, homoclinic 
loop. 

1 Introduction and statement of the main results 



We consider two-dimensional autonomous systems of real differential equations of 
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the form 





where P(x,y) and Q(x,y) are analytic functions defined on an open set U C R 2 . 
Here, the dot denotes, as usual, derivative with respect to the independent variable 
t. The vector field associated to system (JTJ will be denoted X = P(x,y)d x + 
Q(x,y)d y and its divergence is divX := dP/dx + dQ/dy. 

Definition 1 An inverse integrating factor for system ([Tj) in U is a non-locally 
null C 1 solution V : U C M 2 — > M. 0/ £/ie linear partial differential equation 



If V is an inverse integrating factor of system then the zero-set of V, V _1 (0) := 
{(x,y) I V(x, y) = 0}, is composed of trajectories of (jTJ). In fact it is easy to see 
that for any point p G U, if is the orbit of ([I]) that satisfies $(0;p) = p, then 



Thus if V(p) = then p)) = for all t provided that V is defined on $(£; p). 

We consider a certain regular orbit of system (pQ) and we assume the existence 
of an analytic inverse integrating factor in a neighborhood of this regular orbit. 
We show that some qualitative properties of the orbits in a neighborhood of the 
considered regular orbit can be deduced from the known inverse integrating factor. 

We consider a regular orbit <p(t) of system ([1]) and two transversal sections 
£1 and £ 2 based on it. We are going to study the transition map of the flow of 
system ([T]) in a neighborhood of this regular orbit. This transition map is studied 
by means of the Poincare map U : Si — > S2 which is defined as follows. Given 
a point in Si, we consider the orbit of system ([1]) with it as initial point and we 
follow this orbit until it first intersects S2. The map II makes correspond to the 
point in Si the encountered point in S 2 . 

Let (ip(s),ijj(s)) G U, with s G X C R be a parameterization of the regular 
orbit (p(t) between the base points of Si and £2. in particular s can be the time t 
associated to system (pQ). Given a point (x,y) in a sufficiently small neighborhood 
of the orbit (ip(s),i[)(s)), we can always encounter values of the curvilinear coordi- 
nates (s, n) that realize the following change of variables: x(s, n) = <p(s) — ni/j'(s), 
y(s,n) = ip(s) +n(p'(s). We remark that the variable n measures the distance per- 
pendicular to 4>(t) from the point (x, y) and, therefore, n = corresponds to the 
considered regular orbit <f>(t). We can assume, without loss of generality, that the 
transversal section £1 corresponds to £1 := {s = 0} and £2 to £2 := {s = L}, 
for a certain real number L > 0. 

We perform the change to curvilinear coordinates (x, y) 1— > (s, n) in a neigh- 
borhood of the regular orbit n = with s G X = [0,L]. Then, system ([T]) reads 
for 



XV 



VdivX . 



(2) 




(3) 



n = 



N(s,n) 



s = 



S(s,n), 



(4) 



2 



where N(s, 0) = since n = is an orbit and S(s, 0) ^ for s 6 I because 
it is a regular orbit. Therefore, and in order to study the behavior of the orbits 
in a neighborhood of n = 0, we can consider the following ordinary differential 
equation: 

^ = F(s,n). (5) 
as 

We denote by \l/(s; n ) the flow associated to the equation ([5]) with initial condition 
\l/(0;no) = n . In these coordinates, the Poincare map II : Hi — > S 2 between these 
two transversal sections is given by n(n ) = ^(L; n ). 

We assume the existence of an analytic inverse integrating factor V(x, y) in a 
neighborhood of the considered regular orbit 4>(t) of system (OQ). In fact, when 
£i 7^ S 2 and no return is involved, there always exists such an inverse integrating 
factor. It is clear that in a neighborhood of a regular point there is always an 
inverse integrating factor. Applying the characteristics' method, we can extend 
this function following the flow of the system until we find a singular point or a 
return is involved. 

The change to curvilinear coordinates gives us an inverse integrating factor for 
equation ([5]), denoted by V(s,n) and which satisfies 

dV dV dF ~ 

+ ^-F(s,n) = °—V{s,n). (6) 
os on on 

The following theorem gives the relation between the inverse integrating factor 
and the Poincare map defined over the considered regular orbit. 

Theorem 2 We consider a regular orbit (p(t) of system ([I]) which has an inverse 
integrating factor V(x, y) of class C 1 defined in a neighborhood of it and we consider 
the Poincare map associated to the regular orbit between two transversal sections 
II : Si —>■ £2. We perform the change to curvilinear coordinates and we consider 
the ordinary differential equation ([5]) with the inverse integrating factor V(s,n) 
which is obtained from V(x,y). In these coordinates, the transversal sections can 
be taken such that Hi := {s = 0} and E 2 := {s = L}, for a certain real value 
L > 0. We parameterize by a real parameter denoted a. The following identity 
holds. 

V- (£,!%)) = V(0,a)W(a). (7) 



The proof of this result is given in Section [2j 

We remark that if we know an inverse integrating factor for system ([1]), we 
can construct the function V (s,n) and equation ([7]) gives an ordinary differential 
equation for II(cr), which is always of separated variables. Thus, we can determine 
the expression of II (a) in an implicit way and up to quadratures. The first example 
in Section [3] illustrates this remark. 
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We notice that the definition of the Poincare map II is geometric and Theorem 
[2] gives an ordinary differential equation to compute an explicit expression of IT, 
provided that an inverse integrating factor is known. As far as the authors know, 
no such a way of describing the Poincare map has been given in any previous work. 

It is well-known, due to the flow box theorem, that a regular orbit of system 
P]) which is not a separatrix of singular points nor a limit cycle has an associated 
Poincare map conjugated with the identity. There are several results which estab- 
lish that the vanishing set of an inverse integrating factor gives the orbits whose 
associated Poincare map is not conjugated with the identity. In this sense, Giaco- 
mini, Llibre and Viano [13] , showed that any limit cycle 7 C U C IR 2 of system (TTJ) 
satisfies 7 C V -1 (0) provided that the inverse integrating factor V is defined in U. 

Lately it was shown that the zero-set of V often contains the separatrices of 
critical points in U. More precisely, in [3] Berrone and Giacomini proved that, if 
Po is a hyperbolic saddle-point of system ([1]), then any inverse integrating factor 
V defined in a neighborhood of p$ vanishes on all four separatrices of po, provided 
V(po) = 0. We emphasize here that this result does not hold, in general, for non- 
hyperbolic singularities (see [TT] for example). 

We are going to consider regular orbits whose Poincare map is a return map. We 
take profit from the result stated in Theorem [2] in order to study the Poincare map 
associated to a limit cycle or to a homoclinic loop, in terms of the inverse integrating 
factor. Although we have used curvilinear coordinates to state Theorem [2J we do 
not need to use them to establish the results described in the following Theorems 
|3] and [5j The only hypothesis that we need is the existence and the expression of 
an analytic inverse integrating factor in a neighborhood of the limit cycle or the 
homoclinic loop. In particular, we are able to give the cyclicity of a limit cycle 
or of a homoclinic loop from the vanishing multiplicity of the inverse integrating 
factor over the limit cycle or the homoclinic loop. The vanishing multiplicity of an 
analytic inverse integrating factor V(x,y) of system ([1]) over a regular orbit 4>(t) is 
defined as follows. We recall the local change of coordinates x(s, n) = y?(s) — ra//(s), 
y(s,n) = xjj(s) + n(p'(s) defined in a neighborhood of the considered regular orbit 
n = 0. If we have the following Taylor development around n = 0: 

V(x(s,n),y(s,n)) = n m v{s) + 0{n m+1 ), (8) 

where m is an integer with m > 1 and the function v(s) is not identically null, we 
say that V has multiplicity m on (pit). In fact, as we will see in Lemma[6l v(s) 7^ 
for any sGl, and thus, the vanishing multiplicity of V on <f>(t) is well-defined over 
all its points. 
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Let us consider as regular orbit a limit cycle 7 and we use the parameterization 
of 7 in curvilinear coordinates (s,n) with s G [0, L). The Poincare return map IT 
associated to 7 coincides with the previously defined map in which £ := £1 = £2. 
Since the qualitative properties of II do not depend on the chosen transversal 
section £, we take the one given by the curvilinear coordinates and, thus, II(no) = 
no). It is well known that II is analytic in a neighborhood of n Q = 0. We recall 
that the periodic orbit 7 is a limit cycle if, and only if, the Poincare return map II 
is not the identity. If II is the identity, we have that 7 belongs to a period annulus. 
We recall the definition of multiplicity of a limit cycle: 7 is said to be a limit cycle 
of multiplicity 1 if n'(0) 7^ 1 and 7 is said to be a limit cycle of multiplicity m 
with m > 2 if n(n ) = n + f3 m n% + 0(n^ +1 ) with (3 m ^ 0. The following result 
states that given a limit cycle 7 of ([T]) with multiplicity m, then the vanishing 
multiplicity of an analytic inverse integrating factor defined in a neighborhood of 
it, must also be m. 

Theorem 3 Lei 7 be aperiodic orbit of system ([I]) and letV be an analytic inverse 
integrating factor defined in a neighborhood of 7. 

(a) If ^ is a limit cycle of multiplicity m, then V has vanishing multiplicity m 
on 7. 

(b) IfV has vanishing multiplicity m on 7, then 7 is a limit cycle of multiplicity 
m or it belongs to a continuum of periodic orbits. 

The proof of this theorem is given in Section [2j 

The idea to study the multiplicity of a limit cycle by means of the vanishing 
multiplicity of the inverse integrating factor already appears in the work [TJ]. In 
that work the considered limit cycles are semistable and the vanishing multiplicity 
is defined in terms of polar coordinates, which only apply for convex limit cycles. 

Since the Poincare map of a periodic orbit is an analytic function and the 
multiplicity of a limit cycle is a natural number, we obtain the following corollary 
from the previous result. 

Corollary 4 Let 7 be a periodic orbit of system (pQ) and let V be an inverse inte- 
grating factor of class C 1 defined in a neighborhood of 7. We take the change to 
curvilinear coordinates x(s,n) = tp(s) — nip'(s), y(s,n) = ip(s) +rup'(s) defined in 
a neighborhood 0/7. If we have that the leading term in the following development 
around n = 0: 

V(x(s,n),y(s,n)) = n p v(s) + o{n p ) 1 

where v(s) ^0 is such that either p = or p > 1 and p is not a natural number, 
then 7 belongs to a continuum of periodic orbits. 
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We remark that this corollary applies, for instance, in the following case: let us 
consider an invariant curve / = of system ([TJ with an oval 7 such that V/| 7 7^ 
and let us assume that 7 corresponds to a periodic orbit of the system and that 
there exists an inverse integrating factor of the form V = f p g, where g is an nonzero 
function of class C l in a neighborhood of 7 and pel, with either p = or p > 1 
in order to have a C 1 inverse integrating factor in a neighborhood of 7. If we have 
that either p = or p is not a natural number, then we deduce that 7 belongs to 
a period annulus, by Corollary HI On the other hand, applying Theorem [3j in the 
case that p = m G N, we deduce that either 7 belongs to a period annulus or it is 
a limit cycle with multiplicity m. 

A regular orbit <f>(t) = (x(t),y(t)) of ([TJ is called a homoclinic orbit if 4>(t) — > p 
as t — > ±00 for some singular point po. We emphasize that such kind of orbits arise 
in the study of bifurcation phenomena as well as in many applications in several 
sciences. A homoclinic loop is the union T = (j){t) U {po}- We assume that po is a 
hyperbolic saddle, that is, a critical point of system ([TJ such that the eigenvalues 
of the Jacobian matrix DX(p ) are both real, different from zero and of contrary 
sign. We remark that this type of graphics always has associated (maybe only its 
inner or outer neighborhood) a Poincare return map II : E — > S with £ any local 
transversal section through a regular point of T. Moreover, in this work we only 
study compact homoclinic loops that are the a- or cu-limit set of the points in its 
neighborhood, i.e., homoclinic loops with finite singular point and a return map 
different from the identity. This fact implies that T is a compact set, that is, it 
does not have intersection with the equator in the Poincare compactification. In 
fact the same results can be used for some homoclinic loops whose saddle is in 
the equator of the Poincare compactification, but we always assume that the affine 
chart we are working with completely contains the homoclinic loop, that is, the 
homoclinic loops we take into account are compact in the considered affine chart. 

Our goal in this work is to study the cyclicity of the described homoclinic loop 
T in terms of the vanishing multiplicity of an inverse integrating factor. Roughly 
speaking, the cyclicity of T is the maximum number of limit cycles which bifurcate 
from it under a smooth perturbation of ([TJ, see [Tj5] for a precise definition. The 
study of the cyclicity of T has been tackled by the comparison between the Poincare 
return map of the unperturbed system ([1]) and that one of the perturbed system 
in terms of the perturbation parameters. The main result in this study is due to 
Roussarie and appears in [19]. Roussarie presents the asymptotic expression of 
the Poincare map associated to the homoclinic loop T which allows to characterize 
the cyclicity of T. We recall this result in the following Section [2J A good book 
on the subject where these and other results are stated with proofs is [20] • The 
way we contribute to the study of the cyclicity of V is based on the use of inverse 
integrating factors. We take profit from the result of Roussarie to characterize the 
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cyclicity of T as we are able to relate the Poincare return map with the inverse 
integrating factor. 

The Poincare map II : £ — ► S associated to a homoclinic loop V = <fi(t) U{po} is 
given as the composition of the transition map along the regular homoclinic orbit 
4>(t)i which we denote by R, with the transition map in a neighborhood of the 
saddle point po, which we denote by A. We recall that £ is a transversal section 
based on a regular point of T. We have that R is an analytic diffeomorphism 
defined in a neighborhood of the regular part <fi(t) of V and, in this context, it 
coincides with the Poincare map associated along the whole 4>(t) which satisfies 
the identity described in Theorem [2j The map A is the transition map defined 
in a neighborhood of the critical point p^. This map A is characterized by the 
first non-vanishing saddle quantity associated to p . We recall that the first saddle 
quantity is a± = divA'(po) and it classifies the point p between being strong 
(when a>\ ^ 0) or weak (when a% — 0). If po is a weak saddle point, the saddle 
quantities are the obstructions for it to be analytically orbitally linearizable. We 
give the definition of a saddle point to be analytically orbitally linearizable in the 
following Section [21 

In order to define the saddle quantities associated to p , we translate the saddle- 
point po to the origin of coordinates and we make a linear change of variables so that 
its unstable (resp. stable) separatrix has the horizontal (resp. vertical) direction 
at the origin. Let po be a weak hyperbolic saddle point situated at the origin of 
coordinates and whose associated eigenvalues are taken to be ±1 by a rescaling of 
time, if necessary. Then, it is well known, see for instance [T7j, the existence of an 
analytic near-identity change of coordinates that brings the system into: 

fe-i 

x = x + ^2a i x i+1 y i + a k x k+1 y k + ■■■ , 

i=1 (9) 
k-i v ' 

y = -y - ai %iyi+1 ~ bk xkyh+1 + ' ' ' ' 
i=i 

with 0^ — bk 7^ and where the dots denote terms of higher order. The first non- 
vanishing saddle quantity is defined by a k +i := a k — b k . We remark that the first 
non-vanishing saddle quantity can be obtained through an algebraic algorithm. 

We consider a homoclinic loop T = <p(t) U {po} where <fi(t) is a homoclinic orbit 
through the hyperbolic saddle point p . We assume that there exists an analytic 
inverse integrating factor V defined in a neighborhood of the homoclinic loop T. 
The following theorem establishes the relation between the vanishing multiplicity 
of the inverse integrating factor V over the homoclinic orbit <j)(t) and the cyclicity 
of T. We remark that the vanishing multiplicity of the inverse integrating factor 
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also allows to determine the first non-vanishing saddle quantity associated to the 
hyperbolic saddle point p , in case it exists. 

Theorem 5 Let V be a compact homoclinic loop through the hyperbolic saddle po 
of system (CQ) whose Poincare return map is not the identity. Let V be an analytic 
inverse integrating factor defined in a neighborhood ofT with vanishing multiplicity 
m over T. Then, m > 1 and the first possible non-vanishing saddle quantity is 
a m . Moreover, 

(i) the cyclicity ofT is 2m — 1, if a m ^ 0, 

(ii) the cyclicity of T is 2m, otherwise. 

This theorem is proved in Section [2j Theorems [2] and [5] are the main results 
presented in this work. We consider the following situation: we have an analytic 
system (CD) with a compact homoclinic loop T through the hyperbolic saddle point 
Po and whose Poincare return map is not the identity and we aim at determining the 
cyclicity of T. Theorem [5] gives us an algorithmic procedure to solve this problem, 
provided we know an analytic inverse integrating factor V of ([1]) in a neighborhood 
of T with vanishing multiplicity m over it. We compute the m th saddle quantity 
a m at po recalling that all the previous saddle quantities vanish. If a m ^ then 
the cyclicity of T is 2m — 1. Otherwise, when a m = 0, the cyclicity of T is 2m. 

We notice that the determination of the m th saddle quantity can be overcome 
through an algebraic procedure. Thus, our result also applies in the particular 
case in which the saddle of the homoclinic loop is linearizable, that is, the case in 
which a bound for the cyclicity of this graphic cannot be determined through an 
algebraic method. 

We remark that the same result can be applied to a double homoclinic graphic, 
that is, two homoclinic orbits which share the same hyperbolic saddle point po- 
This kind of graphics has also been treated in [T6] . 

This paper is organized as follows. Next section contains the proof of the 
main results, which are stated in this first section. Moreover, in Section [2] several 
additional results related with the problems of studying the cyclicity of a limit cycle 
or a homoclinic loop appear. These results allow the proof of the main results and 
are interesting by themselves as they give the state-of-the-art of the aforementioned 
problems. Much of the results that we present are generalizations of previous ones 
and can be better understood within the context of this work. Last section contains 
several examples which illustrate and complement our results. 
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2 Additional results and proofs 



We consider a regular orbit of system (pQ) and we are going to prove the result 
we have about it, Theorem [21 We consider the aforementioned transition map 
n : Si —>■ £2 along the regular orbit. The key tool to study this Poincare map is 
to change the coordinates of system (CQ) to local coordinates in a neighborhood of 
the regular orbit, that is, to change to the defined curvilinear coordinates (s, n). 
We recall that n = denotes the considered regular orbit and s £ X C R gives 
the transition along the orbits since the two transversal sections can be taken as 
Ei = {s := 0} and E 2 = {s := L}, with L a strictly positive real value. These 
coordinates also allow to give a definition for the Poincare map in terms of the flow 
\I/(s;no) associated to the ordinary differential equation as H(n ) = ^(L; n o)- 
We remark that the expression of the flow \l/(s; no) in a neighborhood of the consid- 
ered orbit no = can be encountered by means of recursive formulae at each order 
of no- This recursive determination of the flow ^(s;n ) allows the study of the 
Poincare map II(too). However, these recursive formulae involve iterated integrals 
which can be, and usually are, very difficult to be computed. The explanation of 
this process and an application to the study of the multiplicity of a limit cycle can 
be encountered in the work |12j . 

We assume the existence of an inverse integrating factor V(x, y) of class C l 
defined in a neighborhood of the regular orbit of system (jTJ). Then, we can con- 
struct an inverse integrating factor for the equation (jSj) which we denote by V(s, n) 
and which satisfies the partial differential equation (jBj). Moreover, the change to 
curvilinear coordinates gives the definition of vanishing multiplicity of V(x, y) over 
the considered regular orbit, as described through the expression (jSJ). 

The following auxiliary lemma gives us that the vanishing multiplicity of V 
over the regular orbit is well-defined on all its points. 

Lemma 6 The function v(s) appearing in (|S]) is different from zero for any s £ X. 

Proof. First, we recall that the relationship between V(x, y) and V(s, n) is V(s, n) = 
V(x(s,n),y(s,n))/(J(s,n)S(s,n)) where J(s,n) is the Jacobian of the change 
(x,y) 1 — > (s,n) to curvilinear coordinates and S(s,n) comes from the time rescal- 
ing t 1 — ► s and is the function defined by system (J3J). We note that J(s,n) = 
(fi'{s) 2 + i^'(s) 2 + 0{n) and thus J(s, 0) 7^ for any s £ X. Moreover, it is clear 
that S(s, 0)^0 for all s £ X. 

The Taylor development of V around n = is of the form 



V(s,n) 



v(s)n m + 0(n' 



m+l 



Here, v(s) = j(s)v(s) where j(s) 



l/( J(s, 0)S(s, 0)) ^ for any s £ X. 
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To end with, we are going to see that v(s) 7^ for any s G X. We recall that, 
by assumption, G X. Since V(s,n) satisfies equation fl§]), we deduce that: 



V (s, n )) = ^ (0, (0; n )) exp | jf — (a, *(«r; n )) da 
We develop this identity around n = and we get that 

u(s)7tf + OK +1 ) = v(0) eocpjjT^foO) dcrj < + 0(n™ +l ) 



(10) 



Since {? (s) ^ 0, we can assume v(0) 7^ 0. Therefore, equating the coefficients of 
in both members of the previous identity, we have the desired result. | 



We have stated all the definitions and notation needed to give the proof of 
Theorem [21 

Proof of Theorem® We have that ^(s;n ) satisfies 



ds 



y s;n 



F (s,^(s;n )) and ^(0;n ) 



no, 



as it is the flow of equation fl5J) with initial condition n . We differentiate the 
previous identities with respect to Uq and we have that: 



d (m 



ds \dn 
Hence, 



dF 
dn 



s,*(s;no)) 



d^ \ 

— - (s;n ) and — (0;n ) 
cm / an 



<9F 

dn 

d^ 



a, ^/(cr; n )) da 



JL ( 99 

da \ dno 



a;n ) 



f|(^o) 



da = In 



dn 



(0-; 



CT=0 



Since ^(0;n ) 
dn 



1, we deduce that 



f dF 

J ~dn ( a '*(°' ;no ^ rf<T = ln 



9* 



Using that the function V(s,n) satisfies (TTU|) . we conclude that 

V(s,*(s;n )) = V-(0,*(0;n )) ^-(s;n ). 

When we take s = L and we recall that the Poincare map is defined as II(no) = 
ty(L;n ), we deduce identity ([7]). | 
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Let us consider as regular orbit a limit cycle 7 and we use the parameterization 
of the whole 7 in curvilinear coordinates with s G [0, L). The value of L can 
be taken the length of the limit cycle 7 and s can be the arc-length parameter 
associated to it. The Poincare return map II associated to 7 coincides with the 
previously defined map in which S := Si = £2. We remark that in curvilinear 
coordinates Si corresponds to the value s = and S 2 to the value s = L and that 
these two sections are equal only in coordinates (x,y). Since 7 is a periodic orbit, 
the change to curvilinear coordinates and all the analytic functions in (x, y) in a 
neighborhood of 7 need to be L-periodic in s. 

We use that V(s, n) is L-periodic in s and we develop relation (JTj) in a neigh- 
borhood of a = in order to get the following result. We recall that the point 
cr = corresponds to the periodic orbit 7 and that it is, therefore, a fixed point of 
its Poincare map. 

Proposition 7 Let us consider the following Taylor development of the analytic 
Poincare map 11(a) associated to a periodic orbit 7, around a = 

( e^a + 0(a 2 ) with ft^O, 

n(o-) = I a + (3 k a k + 0(a k+1 ) with (3 k ^ 0, k > 1 , 

^ cr otherwise. 

Let m > be the vanishing multiplicity of an analytic inverse integrating factor V 
on 7. 

(a) lfm = 0, then (3 k = for all k > 1. 

(b) If m > 1, then either m is the lowest subindex such that (3 m 7^ or (3k = 
for all k > 1 . 

Proof. Let us recall the Taylor development of the function V(s,no) in a neigh- 
borhood of n Q = 0: V"(s,no) = J2i>o^i( s ) n o- Since the function V(x, y) is not 
locally null, we have that there exists an integer m with m > which is the lowest 
index with v m (s) 7^ for any s 6l from Lemma [6j By definition, this value of m 
is the vanishing multiplicity of V(x,y) on the considered regular orbit. 

Since V(x, y) changed to local coordinates is L-periodic in s, we have that the 
function V(s, n) is also L-periodic in s. Then, we can write V (0, a) = J2i>o v i a% 
and\/(L,n(cr)) = Ei^o^nW*- 

If (3x 7^ 0, we evaluate identity (j7]) in cr = and we get that (L, 0) = 
e^V (0, 0) which implies that vq = and, thus, the vanishing multiplicity of V(x, y) 
on the regular orbit is at least 1. We develop identity © in powers of cr and the 
lowest order terms in both members of the equality correspond to a m . The equation 
for their coefficients is: z/ m e m/31 = ^ m e ft . Since (3% 7^ 0, we deduce that m = 1. 
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Let k > 1 be the lowest subindex such that fa 7^ 0. We have that n(cr) = 
a + faa k + 0(a k+l ). We subtract V (0, a) from both members of ([7D and we get 
the following relation: 

(U(aY - a 4 ) = U + X>i<r') (kfacr^ 1 + 0(a fc )) . (11) 

The left hand side of (TTTT) has order at least o~ k which implies that uq = 0, in order 
to have the same order in both members. Therefore, m > 0. We have that the 
lowest order terms in both sides of (TTTT) correspond to C r m+fc_1 and the equation of 
their coefficients is: mfav m = k fa v m , which implies that k = m. | 

Proof of Theorem® It is a straightforward consequence of the previous Propo- 
sition [3 We have preferred to state this proposition as we are going to use the 
same reasoning for the transition map associated to a homoclinic orbit. | 

Let us now consider a homoclinic loop T = <f>(t) U {p } whose critical point p 
is a hyperbolic saddle-point. We will denote by A and fi the eigenvalues associated 
to the Jacobian matrix DX(p ) with the convention fi < < A. We associate to 
Po its hyperbolicity ratio r = — fi/\. We say that the singular point p is strong if 
divX(p ) 7^ (equivalently r^l) and it is weak otherwise. The hyperbolic saddle 
Po is called p : q resonant if r = q/p G Q + with p and q natural and coprime 
numbers. 

A homoclinic loop T is called stable (unstable) if all the trajectories in some 
inner or outer neighborhood of T approach r as t — > +00 (t — > —00). In the 
investigation of the stability of a homoclinic loop T of system ([TJ through a saddle 
Po the quantity 

ai = divX( Po ) (12) 

plays an important role. In short, it is well known, see for instance p. 304 of 
[TJ, that T is stable (unstable) if a\ < (a± > 0). For this reason such kind of 
homoclinic loops T are said to be simple if ot\ 7^ and multiple otherwise. 

The cyclicity of V is linked with its stability. It is notable to observe that, 
in the simple case a± 7^ 0, the stability of T is only determined by the nature of 
the saddle point itself. Andronov et al. [TJ proved that if a% 7^ 0, the possible 
limit cycle that bifurcates from T after perturbing the system by a multiparameter 
family has the same type of stability than T, hence this limit cycle is unique for 
small values of the parameters. After that, Cherkas [8] showed that, if ati — 
and the associated Poincare map of the loop T is hyperbolic, then the maximum 
number of limit cycles that can appear near T perturbing the system in the C l class 
is 2. We recall that a real map is said to be hyperbolic at a point if its derivative at 
the point has modulus different from 1. In [19] Roussarie presents a generalization 



12 



of these results which determine the asymptotic expression of the Poincare map 
associated to the loop T. 

The Poincare map II associated to T is defined over a transversal section £ 
whose base point is a regular point of T. We parameterize the transversal section 
£ by the real local coordinate a. The value a = corresponds to the intersection 
of £ with T and o > is the side of T where n(cr) is defined. Roussarie's result is 
twofold: on one hand the asymptotic expansion of n(<r) is determined and, on the 
other hand, the cyclicity of T is deduced from it. 



Theorem 8 | jl Qj | Let us consider any smooth perturbation of system ([I]) . Then we 
have: 

(i) If r 7^ 1 (equivalently, «i ^ 0), then H(a) = ca r {l + o(l)) with c > and at 
most 1 limit cycle can bifurcate from T. 

(ii) If Oi\ = and Pi ^ 0, then IT((x) = Gr x a + o(er) and at most 2 limit cycles 
can bifurcate from V . 

(Hi) If cti = Pi = for % = 1,2, ... ,k with k > 1 and a^+i ^ 0, then n(cr) = 
a + ak+ic k+1 log cr + o(cr fc+1 loger) and at most 2k + 1 limit cycles can bifurcate 
from T. 

(iv) If cti = Pi = oik = for i = 1, 2, . . . , k — 1 with k > 1 and Pk ^ ; then 
n(cr) = o + Pk& k + o(a k ) and at most 2k limit cycles can bifurcate from V. 

(v) If cti — Pi — for all % > 1, then Tl(o~) = o and the number of limit cycles 
that can bifurcate from T has no upper bound. 

The values «j, with i > 1, are the saddle quantities associated to po and are 
evaluated by a local computation. The values Pi, with i > 1, are called separatrix 
quantities of V and correspond to a global computation. The determination of the 
«i can be explicitly done through the algebraic method of normal form theory near 
Pq. Additionally, Pi = j T divX dt, and a more complicated expression for p 2 can 
be encountered in [16] and it involves several iterated integrals. On the contrary, 
as far as we know, there is no closed form expression to get Pi for i > 3. We remark 
that our result provides a way to determine the cyclicity of the homoclinic loop T 
without computing any 

A graphic f = U k =l (fii(t) U {pi, . . . ,pk} is formed by k singular points pi, ■ ■ ■ ,Pk, 
Pk+i = Pi and k oriented regular orbits 4>i(t), . . . , 4>k(t), connecting them such that 
4>i(t) is an unstable characteristic orbit of Pi and a stable characteristic orbit of 
Pi+\. A graphic may or may not have associated a Poincare return map. In case it 
has one, it is called a polycycle. Of course, the homoclinic loops that we consider 
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are poly cycles with just one singular point. In [3J, the case of a system §Q with a 
compact polycycle T whose associated Poincare map is not the identity and whose 
critical points pi,Pz, ■ ■ ■ ,Pk are all non-degenerate is studied and it is shown that 
if there exists an inverse integrating factor V defined in a neighborhood of F, then 
r C y -1 (0). A strong generalization of this result is given in [11] , where it is showed 
that any compact polycycle contained in IA with non-identity Poincare return map 
is contained into the zero-set of V under mild conditions. More concretely, Garcia 
and Shafer give the following result. 



Theorem 9 jllj Assume the existence of an analytic inverse integrating factor 



V defined in a neighborhood M of any compact polycycle V of system ([[]) with 
associated Poincare return map different from the identity. Then, T C V _1 (0). 

In this result is also given for an inverse integrating factor V with lower 
regularity than analytic and assuming several conditions. In fact, a consequence 
of their results is that if we consider a compact homoclinic loop T whose Poincare 
map is not the identity and such that there exists an inverse integrating factor V 
of class C l defined in a neighborhood of T, then T C V -1 (0). 

Along this paper, we will work with inverse integrating factors V(x, y) analytic 
in a neighborhood of a compact homoclinic loop T = <fi(t) U {p^} with associated 
Poincare map different from the identity. Hence, T C V -1 (0). 

Regarding the existence problem, it is well known that the partial differential 
equation (|2J) has a solution in a neighborhood of any regular point, but not neces- 
sarily elsewhere. In [6] it is shown that if system ([1]) is analytic in a neighborhood 
M C U of a critical point that is either a strong focus, a nonresonant hyperbolic 
node, or a Siegel hyperbolic saddle, then there exists a unique analytic inverse 
integrating factor on M up to a multiplicative constant. Hence, we only have at 
the moment the local existence of inverse integrating factors in a neighborhood 
of convenient singular points. But, we shall need to know a priori whether there 
exists an inverse integrating factor V defined on a neighborhood of the whole ho- 
moclinic loop T of system ([I]). This is a hard nonlocal problem of existence of 
global solutions of the partial differential equation ([2]) for which we do not know 
its answer. We describe one obstruction to the existence of an analytic inverse 
integrating factor defined in a neighborhood of certain homoclinic loops. 

Proposition 10 Suppose that system (CQ) has a homoclinic loop T through the 
hyperbolic saddle point po which is not orbitally linearizable, p : q resonant and 
strong (p^ q). Then, there is no analytic inverse integrating factor V(x,y) defined 
in a neighborhood ofT. 

Proof. Let us assume that there exists an analytic inverse integrating factor V(x, y) 
defined in a neighborhood of Y. Let us consider f\(x,y) = and f^x.y) = the 
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local analytic expressions of each of the separatrices associated to po, where the 
subindex denotes the corresponding eigenvalue. As we will see in Proposition [T2] 
(Hi), we have that V(x, y) factorizes, as analytic function in a neighborhood of po, 
in the form V = fl +kq fa +kPu with u(po) ^ and an integer k > 0. Since V(x, y) 
is defined on the whole loop T, it needs to have the same vanishing multiplicity 
on each separatrix at po, which implies that k = 0. Thus, we have that the local 
factorization of V in a neighborhood of p is V = fx / M u. 

Let us take local coordinates in a neighborhood of po, which we assume to be 
at the origin. We recall, see [2] and the references therein, that given an analytic 
system x — Xx + ■ • -, y = \iy + ■ • • near the origin with = —q/p G Q~ with 
p and q natural and coprime numbers (p : q resonant saddle) then the system has 
two analytic invariant curves passing through the origin. Moreover, it is formally 
orbitally equivalent to 

X = pX [1 + 5(U e + aU 21 )] , Y = -qY , (13) 

with U = X q Y p , a e R, I an integer such that £ > 1 and 5 G {0, ±1}. The normal 
form theory ensures that po is orbitally linearizable if, and only if, 5 = 0. We are 
assuming that po is not orbitally linearizable and, thus, 5^0. 

Let us consider V(X, Y) the formal inverse integrating factor of system ffl3|) 
constructed with the transformation of V(x, y) with the near-identity normalizing 
change of variables divided by the jacobian of the change. We have that V(X, Y) = 
XYu(X, Y) where u is a formal series such that w(0,0) 7^ 0. We can assume, 
without loss of generality, that u(0, 0) = 1. Easy computations show that if u is 
constant, we have that V(X, Y) cannot be an inverse integrating factor of system 
(1T3"|) with 5 7^ 0. If we have that u is not a constant, we develop uasa formal series 
in X and Y and we can write u = 1 + V S (X, Y) + ■ ■ ■ where the dots correspond to 
terms of order strictly greater than s in X and Y and V S (X,Y) is a homogenous 
polynomial in X and Y of degree s. We consider the partial differential equation 
satisfied by V: 

pX[l + 6(Ut + aU*)]^- q Y§ = 

= (p - q + p6(l + lq) U e + a5p(l + 2£q) U n ) V. 
We equate terms of the same lowest degree, we deduce that s = £(p + q) and that: 

pX m ~ qY w = b ^ xlqylv - 

Easy computations show that the general solution of this partial differential equa- 
tion is V S (X,Y) = 5£ 2 pq 2 U e In X + G{U) where G is an arbitrary function. We 
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deduce that this partial differential equation has no polynomial solution V S (X,Y) 
unless 5 = 0. 

We conclude that the existence of such an inverse integrating factor V(X, Y) 
implies that 5 = in contradiction with our hypothesis. | 

We observe that homoclinic loops considered in Proposition [TU1 are simple, since 
the hyperbolic saddle point is strong. As we have already mentioned in the intro- 
duction, it is well-known that its cyclicity is 1. Therefore, the obstruction to the 
existence of an analytic inverse integrating factor in a neighborhood of these ho- 
moclinic loops is not relevant in the context of the bifurcation theory. 

The following example provides a realization of the thesis stated in Proposition 
[TOj that is, we illustrate the existence of homoclinic loops where no analytic inverse 
integrating factor can exist in a neighborhood of it. We consider the system 

x = —x + 2y + x 2 , y = 2x — y — 3x 2 + —xy , (14) 

studied in [I]. The origin is a strong saddle because it has eigenvalues \i = — 3 
and A = 1 and, hence r = 3 ^ 1. Moreover, system (|14p possesses a homoclinic 
loop T through the origin contained in the invariant algebraic curve f(x,y) = 
x 2 (l — x) — y 2 = and having inside a focus. After a linear change of variables, 
we write system (1141) with Jordan linear part as 



o 2 i o ^2 • ^2, 1 2 

x = —6x — —x + 2xy — —y , y = y — -x + xy — -y . 

A computation using the normal form method shows that the above system is 
conjugated to the system 

86579 4 

x = —6x , y = y xy , 

y y 248832 y 

up to homogeneous degree 15 included. Therefore we conclude that system (fT4j) is 
not formally orbitally linearizable near the origin. Thus, applying Proposition (TDl 
it exists no analytic inverse integrating factor V(x,y) for system (THl) defined in a 
neighborhood of V. 

We have used the method described in [7] to show that system fTl4|) has no 
irreducible invariant algebraic curve except x 2 (l — x) — y 2 = 0, nor exponential 
factors. Therefore, from the Darboux theory of integrability it, can be shown that 
this system has no Liouvillian first integral. For the definitions and results related 
to invariant algebraic curves, exponential factors and the Darboux theory of inte- 
grability, see [7] and the references therein. 
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In this paper we will always assume as a hypothesis the existence of an analytic 
inverse integrating factor V defined on a neighborhood of T. Under this condition, 
we now refer to the uniqueness problem. In the last section of this work, we present 
several examples of differential systems with a homoclinic loop which satisfy all our 
hypothesis and in which an explicit expression of an analytic inverse integrating 
factor is given. 

Proposition 11 Consider the analytic system having a compact loop F with a 
Poincare return map different from the identity. Then, assuming the existence of 
an analytic inverse integrating factor V(x, y) of ([T]) defined in a neighborhood of 
T, we have that V(x,y) is unique up to a multiplicative constant. 

Proof. Assume that we have two different analytic inverse integrating factors V and 
V defined in a neighborhood M of V. From Theorem O it follows that V C V _1 (0) 
and T C U _1 (0). Let m and fh be the multiplicities of vanishing of V and V on 
T, respectively. We can assume m > m and we can construct for system ([1]) the 
first integral H(x,y) = V/V which is analytic in M. Lemma [6] ensures that this 
quotient has no poles in N ' . The existence of H is in contradiction with the fact 
that T is an a- or Lj-limit set with a return map different from the identity. | 

As a corollary of the proof of this proposition we have that if there exists an 
analytic inverse integrating factor defined in a neighborhood of a limit cycle, then 
it is unique up to a multiplicative constant. 

An easy reasoning shows that if we assume the existence of an analytic inverse 
integrating factor V(x, y) defined in a neighborhood of a compact loop T, through 
the singular point po and whose Poincare return map is not the identity, and if 
there exists a unique formal inverse integrating factor V(x, y) in a neighborhood 
D of po, then V(x, y) = V(x,y) for all (x, y) G D up to a multiplicative constant. 

Let n(a) be the Poincare return map associated to a homoclinic loop V and 
defined over a transversal section S. Then, II is the composition II = R o A where 
R is an analytic diffeomorphism defined in a neighborhood of the regular part 4>(t) 
of T and A is the transition map defined in a neighborhood of the critical point 
Pq. We parameterize E by a real parameter a > and o = corresponds to the 
base point. In our notation a denotes both a point in E and the corresponding 
real value which parameterizes it. Theorem [2] relates the inverse integrating factor 
and the diffeomorphism defined in neighborhood of any regular orbit by the flow of 
system ([1]). Let us consider a homoclinic orbit as regular orbit and we remark that 
Proposition [7| also applies in the following sense. In case we consider a homoclinic 
orbit, the limiting case in which the two sections tend to the saddle po with Hi 
following the unstable separatrix and £2 the stable separatrix, gives that this 
Poincare map coincides with the aforementioned regular map R. We remark that, 
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if X = (a, b) covers the whole homoclinic orbit in the sense of the flow, this limiting 
case corresponds to go from £ 1; based on a value of s such that s -> a + , to E 2 
based on a value of s such that s — > b~. 

We note that, since V(x, y) is a well-defined function in a neighborhood of 
the whole T, the function V (s, n) takes the same value in the limiting case, i.e., 
when s tends to the boundaries of X. We define V (a, n) = lim s ^ a + V (s, n) and 
V (b, n) = lim s ^ b - V (s, n) and these limits exist from the same reasoning. More- 
over, V(a,n) = V (b,n). Thus, the result stated in Proposition [7] is also valid to 
study the regular map R associated to a homoclinic orbit defined with sections in 
the limiting case. 

In order to relate inverse integrating factors and the Poincare return map as- 
sociated to a homoclinic loop, we need to study the local behavior of the solutions 
in a neighborhood of the critical saddle point pq. 

The following proposition establishes some relationships between the vanishing 
multiplicity of V on the separatrices of a hyperbolic saddle p and the nature of 
Po, provided that V(p ) = 0. 

Proposition 12 Let V(x, y) be an analytic inverse integrating factor defined in a 
neighborhood of a hyperbolic saddle po with eigenvalues \i < < A of an analytic 
system (pQ). Let us consider f\(x,y) = and f^(x,y) = the local analytic expres- 
sion of each of the separatrices associated to po, where the subindex denotes the 
corresponding eigenvalue. Then the next statements hold: 

(i) If Pq is strong, then V(p ) = 0. 

(ii) If Po is nonresonant, then V = f\ f^u with u(po) ^ 0. 

(Hi) If Pq is p : q resonant and strong, then V = fl +kq f]^ kv u with u(p ) ^ and 
an integer k > 0. 

(iv) If po is weak and V(po) = ; then V = f™ f™ u with u(po) ^ and m a 
natural number with m > 1. 

Proof. Statement (i) is clear from the definition (j2j) of inverse integrating factor. 

The other three statements are based on a result due to Seidenberg [2T], see 
also [7]. Since V is analytic, we can locally factorize near po as V = f™ 1 f™ 2 u with 
u(po) 7^ and m« nonnegative integers. Then, we have divA'(po) = miX + m2/i. 
Hence, (mi — 1)A + (m.2 — = and statements (ii)-(iv) easily follow. | 

Let us consider an analytic system (CQ) where po is a weak hyperbolic saddle. 
By an affine change of coordinates and rescaling the time, if necessary, the system 
can be written as 

x = x + f(x,y),y = -y + g(x,y), (15) 
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where / and g are analytic in a neighborhood of the origin with lowest terms at 
least of second order. It is well known the existence of a formal near-identity 
change of coordinates (x, y) i— > (X, Y) — (x + ■ ■ ■ , y + ■ ■ •) that brings system ffl5|) 
into the Poincare normal form 



X = X 



i>l 



,Y = -Y 



1 + 



i>l 



(16) 



From this expression, we see that system (Tl5|) has an analytic first integral in a 
neighborhood of the saddle if and only if the saddle quantities oti + \ := — hi are 
zero for all z > 1, see [T7] and references therein for a review. In particular we 
observe that system (TTBT) has a 1:1 resonant saddle at the origin. In [15], it is proved 
that a planar dynamical system is analytically orbitally linearizable at a resonant 
hyperbolic saddle, that is whose hyperbolicity ratio r is a rational number, if and 
only if it has an analytic first integral in a neighborhood of the saddle. We recall 
that a saddle is analytically orbitally linearizable if there exists an analytic near- 
identity change of coordinates transforming the system to a local normal form such 
that Y/X — —rY/X in a neighborhood of the saddle. 

We say that system (Tl5l) is analytically orbitally linearizable (integrable) at the 
origin if there exists an analytic near-identity change of coordinates transforming 
the system to 

X = Xh{X, Y) , Y = -Yh(X, Y) , 

with h(0, 0) = 1. In this case, V(X, Y) = X k Y k h(X, Y) is a 1-parameter family 
of analytic inverse integrating factors for any integer k > 1. 

In [I], Brjuno shows that any resonant hyperbolic saddle point of an analytic 
system is analytically orbitally linearizable if and only if it is formally orbitally 
linearizable. In particular, this fact means that either there exists at least one 
saddle quantity aj different from zero or the system becomes analytically orbitally 
linearizable. Moreover, we can prove the following result. 

Theorem 13 Let us consider an analytic system ([I]) with a hyperbolic weak saddle 
p whose separatrices are locally given by fi(x,y) = and f-i(x, y) = 0. Let 
V = f™ fH\ u be an analytic inverse integrating factor defined in a neighborhood 
of po such that u(p ) ^ and m is a nonnegative integer. If p is not analytically 
orbitally linearizable, then m > 2 and the first nonvanishing saddle quantity is a m . 

Proof. We can always assume that system ([T|) takes the form ([15]) where the 
hyperbolic weak saddle po is at the origin and the expression of the separatrices 
takes the form fi(x,y) = x + o(x,y), f^\(x,y) = y + o(x,y). As a result of 
Proposition [T21 the analytic inverse integrating factor reads for V(x, y) = f™f™iU 
with nonnegative integer m > and w(0, 0) ^ 0. 
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If m = 0, then V(po) 7^ and there exists an analytic first integral defined on 
a neighborhood of po for system ([T]). Hence, po is analytically orbitally linearizable 
and the saddle quantities of po satisfy = for any i > 1. 

If m = 1, then V = fif-iu with w(po) 7^ 0. Then, by statement (iii) of Theorem 
5.10 of [9], we conclude that p Q is also analytically orbitally linearizable and the 
saddle quantities of p satisfy aij = for any i > 1. 

If m > 2, we are going to prove that the saddle quantities at the origin satisfy 
«i = for i = 1, 2, . . . , m — 1 and that a m 7^ 0. We do a formal near-identity 
change of coordinates (x, y) 1— > (X, K) transforming system (|T5|) into the Poincare 
normal form (fT6l) . In these normalizing coordinates, system (fT6j) has the formal 
inverse integrating factor 

V(X,Y)=XYY J ^ + i(XYY • (17) 

We observe that functions which define system (fT6l) give V"(X, F) = X/X — Y/Y. 
We remark that V"(X, Y) is a univariate function of the variable XY . Since there is 
at least one nonzero saddle quantity, system (fT6l) has no formal first integral near 
the origin. In particular, f[T7|) is the unique (up to multiplicative constants) formal 
inverse integrating factor of system ffTBl) . We note that performing the above for- 
mal near-identity change of coordinates to the inverse integrating factor V(x, y) we 
get a formal inverse integrating factor of system (TL6]) of the form X m Y m u(X, Y), 
with m(0, 0) 7^ 0. Therefore, from uniqueness, we must have that it coincides with 
the expression (1T7I) . Thus, we have «j = for i = 1, 2, . . . , m—1 and that a m 7^ 0. | 

We have analyzed the behavior of the flow near a hyperbolic saddle point po 
when an analytic inverse integrating factor V is known. In particular, we shall see 
that the transition map A near po can be studied by the vanishing multiplicity of 
V in the separatrices of Pq. Let us consider two transversal sections E 2 and Si in 
a sufficiently small neighborhood of po, where S 2 (resp. X^) is based on a point 
over the stable (resp. unstable) separatrix of po. We parameterize £2 by a real 
parameter a > and a = corresponds to the base point. We recall that the 
transition map A : £2 —* £1 is defined as A(cr) = $(r(<r);cr) D Si where $(t;cr) 
is the flow associated to system ([T|) with initial condition the point a G S 2 and 
r((x) is the minimal positive time such that the trajectory through a intersects S x . 
The explicit asymptotic expansion of A(er) was given by Dulac [TU] in terms of the 
hyperbolicity ratio r and saddle quantities of po as follows: 

( ca r (l + o(l)) ,c>0 if r 7^ 1 , 
= < cr + ttfcO" fc lncr + o(cr fc lncr) if a\ = ■ ■ ■ = otk-i = 0, at 7^ , 
^ a if a,i — 0, for all i > 1. 

Since the Poincare return map II(cr) associated to the compact homoclinic loop 
T is the composition II = R o A and is different from the identity, we only have 
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the four possibilities (i)-(iv) described in Theorem [HI 

In the following theorem, the cyclicity of Y denotes the maximum number of 
limit cycles that bifurcate from T under smooth perturbations of ([I]). 

Theorem 14 Let Y be a compact homoclinic loop through the hyperbolic saddle po 
of system §Q whose Poincare return map is not the identity. Let V be an analytic 
inverse integrating factor defined in a neighborhood ofY with vanishing multiplicity 
m over Y . Then the following statements hold: 

(a) m > 1. 

(b) If Pq is strong, then m = 1 and the cyclicity of Y is I. 

(c) If pq is weak, then: 

(c.l) If Pq is not analytically orbitally linearizable, then m > 2, cti = fa = 
for i = 1,2, . . .,m — 1 and a m ^ 0. Moreover, the cyclicity of Y is 
2m - 1. 

(c.2) Ifp is analytically orbitally linearizable, then fa — fa = • • • = /3 m -i — 
and f3 m ^ 0. Moreover, the cyclicity ofY is 2m. 

Proof. Statement (a) is a straight consequence of the fact that Y C y _1 (0), see 
Theorem 

Statement (b) follows taking into account statements (ii) and (iii) of Proposition 
[T2l where we recall that the vanishing multiplicity of V in each of the separatrices 
of po needs to be the same since T is a loop. In particular, we observe that the 
value of k in statement (iii) of Proposition [121 must be zero by the same argument. 

We assume that p is weak. We may have that all its associated saddle quan- 
tities are zero or that there is at least one saddle quantity different from zero. If 
this last case applies, the first non-vanishing saddle quantity result of 

Theorem [13j Therefore, we compute the m th saddle quantity associated to po, 
knowing than the previous saddle quantities need to be zero and we can determine 
if po is analytically orbitally linearizable (if a m = 0) or not (if a m ^ 0). The case 
(c.l) in the theorem corresponds to a m ^ and the case (c.2) to a m = 0. 

The case (c.l) is a consequence of Theorem [13] and part (b) of Proposition 
[7[ More concretely, since po is not analytically orbitally linearizable, then using 
Theorem [TBI we have m > 2, = for % = 1, 2, . . . , m — 1 and a m ^ 0. In addition, 
by part (b) of Proposition [7J we get fa = for i = 1, 2, . . . , m— 1. We remark that, 
in this case, either (3 m ^ or fa = for all k > 1. In any case, from statement 
(iii) of Theorem [HI the cyclicity of Y is 2m — 1. 

Finally, the proof of (c.2) works as follows. Since po is analytically orbitally 
linearizable, then at = for all k > 1. By hypothesis, the Poincare map is not the 
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identity, so (3k = for all k > 1 is not possible. Hence, by part (b) of Proposition 
[7J we get that /3j = for % = 1, 2, . . . , m — 1 and /3 m ^ 0. Thus, using statement 
(iv) of Theorem [HI the cyclicity of V is 2m. | 

Proof of Theorem^ The thesis of Theorem [5] is a corollary of Theorem [TH | 



3 Examples 

This section contains several examples which illustrate and complete our results. 
The first example shows how to give an implicit expression of the Poincare map 
associated to a regular orbit via identity The second example is given to show 
the existence of analytic planar differential systems with a compact homoclinic 
loop T whose cyclicity can be determined by the vanishing multiplicity of an in- 
verse integrating factor defined in a neighborhood of it. The third example targets 
to show that there is no upper bound for the number of limit cycles which bifurcate 
from a compact homoclinic loop whose Poincare return map is the identity. The 
system considered in this example has a numerable family of inverse integrating 
factors such that for any natural number n there exists an analytic inverse inte- 
grating factor whose vanishing multiplicity on the compact homoclinic loop is n. 

Example 1. The following system is studied in |5] where a complete description 
of its phase portrait in terms of the parameters is given using the inverse integrating 
factor as a key tool. The system 

x = Xx — y + XrriiX 3 + (m 2 — mi + mim 2 )x 2 y + Xm 1 m 2 xy 2 + m 2 y 3 , 
y = x + Xy — x 3 + Xm\x 2 y + (mim 2 — mi — l)xy 2 + Amim 2 ?/ 3 , 

where A, mi and m 2 are arbitrary real parameters, has the following inverse inte- 
grating factor 

V(x,y) = (x 2 + y 2 ) (1 + mil 2 + mim 2 y 2 ). (19) 

In [5] it is shown that if A ^ 0, mi < 0, mi ^ — 1 and m 2 > 0, then the ellipse 
defined by H-m 1 x 2 + m 1 m 2 ?/ 2 = 0, and which we denote by 7, is a hyperbolic limit 
cycle of system ( Fl8|) . Moreover, with the described values of the parameters, we 
have that the origin of coordinates is a strong focus whose boundary of the focal 
region is 7. Moreover, since the vanishing set of V is only the focus at the origin 
and the limit cycle 7, we deduce by Theorem [HJ that the region, in the outside of 
7 in which the orbits spiral towards or backwards to it, is unbounded. 

We remark that we reencounter that this limit cycle is hyperbolic using that the 
vanishing multiplicity of the inverse integrating factor over it is 1. In this example 
we are not concerned with the multiplicity of the limit cycle but with the Poincare 
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return map associated to it. We are going to use the ordinary differential equation 
stated in (|7|) to give an implicit expression for the Poincare map associated to 7 in 
system ( ITS]) . We can parameterize the ellipse by (a/tT?^ cos(s), sin(s))/^/-m 1 m 2 
with s € [0, 27r) and we perform the corresponding change to curvilinear coor- 
dinates: (x,y) 1— > (s, n) with 2 = (y 7 ^ — n) cos(s)/y/— m\m2 and y = (1 — 
wfr^n) sin(s) / ^/—mxm^- We obtain an ordinary differential equation of the form 
(JH]) which describes the behavior of the solutions of system (ITS]) in a neighborhood 
of the limit cycle (for n = 0). The corresponding inverse integrating factor V(s,n) 
is a 27r-periodic function in s and it satisfies that: 

mi n { n 



rr, n \ mi n (n - V m 2j I" - <V m 2 

1/(0, n) = — ^ 7= 

n 1 — 2riy/m2 + m.2 + mirri2 



The ordinary differential equation for this system reads for V (2ir, n(er)) = 
(0, a) n'(cr), where cr can be taken as a real parameter of the section X := 
{s = 0} and such that a = corresponds to the limit cycle. We observe that the 
point cr = v /m2 corresponds to the focus point at the origin of system (TT8]) . This 
differential equation can be written in the following Pfaffian form: 



da dU 



V(0,a) V{2tt,U) 



whose integration gives: 



(n(q)(n(a)-2 v ^)) (1+mi)/2 = (a(a-2 V m- 2 )f +m ^ 2 

(n((j) - v/m^)" 11 ° (cr - v^)" 11 ' ^ ' 

where fco is an arbitrary constant of integration. We remark that if ko — 1, then 
the function n(cr) := cr satisfies the implicit identity ( [20] ). 

We observe that the integration of (j7]) always gives rise to a Pfaffian form of 
separated variables. Moreover, in the case that a return is involved we have that 
this Pfaffian form is symmetric in IT and cr, because the function V(s, n) needs to 
be L-periodic in s. 

We are interested in the Poincare map associated to the limit cycle 7. Since 
it is a hyperbolic limit cycle, we have that II(er) = e^cr + 0(a 2 ) with (3\ ^ 0, 
as described in Proposition [7] Using this expression and identity ([20]) . we deduce 
that k = exp (j3i (mi + l)/2). Moreover, it can be shown that /3i = —2\T, where 
T is the minimal positive period of 7 and takes the value T = 2nmi/(l + mi) if 
mi < —1 and T = — 27rmi/(l + mi) if —1 < mi < 0. Therefore, ko = e~ 2Xmin if 
mi < — 1 and ko = e 2Ami7r if —1 < vii\ < 0. We note that the inverse integrating 
factor determines the multiplicity of the limit cycle 7 but not its stability. 
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Example 2. We fix an integer number m with m > 1 and we consider the 
algebraic curve / = with f(x, y) = y 2 — (1 — x) 2 (l + x). We have that / = has 
an oval in the range — 1 < x < 1 with a double point in (1,0). We denote by V 
this oval. Let us consider the planar differential system: 



x 



( [(1 - m)g(x, y) + f(x, y)^ 1 ] ^ + f(x, y)^j (x 2 + y 2 )q(x, y) - 
f{x, y) m \2{x + y)q(x, y) + (x 2 + y 2 )- 



dy 

[(1 - m)g(x, y) + f(x, y)^ 1 ] ^ + f(x, y)^\ (x 2 + y 2 )q(x, y) + 
f(x,y) m (2(x-y)q(x,y) + (x 2 + y 2 ) 



dx 

(21) 

where g(x,y) and q(x,y) are polynomials such that the algebraic curve g = 
does not intersect the oval T and the algebraic curve q = does not contain any 
point in the closed region bounded by T. This system has T as homoclinic loop 
where the critical point (1, 0) is a hyperbolic weak saddle. Moreover, the origin 
of this system is a strong focus. It can be shown that the function V(x, y) = 
(x 2 + y 2 )f(x,y) m q(x,y) is an inverse integrating factor, which is analytic in the 
whole affine plane. Since the origin is a strong focus and V is analytic in R 2 , we 
have by Theorem [9], that the boundary of the focal region must be contained in the 
zero-set of V. We conclude that the boundary of this focal region needs to be T 
and, therefore, we have that T is a compact homoclinic loop through a hyperbolic 
saddle and whose associated Poincare return map is not the identity. Moreover, 
we deduce that the vanishing multiplicity of V over T is m. The following function 

H{x,y) = (x + iyY-'ix -iy) 1+i f(x,y)e 9{x ' y)/f ^ yr ~ 1 q(x,y) 

is a first integral of system (|2ip . In the work [IB] it is shown that this form of 
a first integral implies that the inverse integrating factor is polynomial and, thus, 
well-defined over all the real plane. When m = 1, H(x,y) provides an analytic 
first integral defined in a neighborhood of the saddle point (1, 0) and we deduce 
that the saddle quantities associated to the critical point (1,0) are all zero. Thus, 
using Theorem [T3J we conclude that when m = 1 the cyclicity of V in system (12 ip 
is 2. 

When m > 1, using Theorem [TU we deduce that the cyclicity of F in system 
f l2T|) is 2m — 1 if a m ^ and 2m if a m = 0. 

Example 3. Let us consider the following planar differential system 

x = -2y, y = -2x + 3x 2 , (22) 
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which has a homoclinic loop T contained in the invariant algebraic curve / = 
with f(x, y) := y 2 — x 2 + x 3 . Since the origin (0, 0) is the saddle contained in T and 
it is a hyperbolic saddle, we have that T has associated a Poincare return map, 
which is the identity as the system is Hamiltonian with H = f. In addition to the 
hyperbolic saddle at the origin, the system possesses a singular point of center type 
in the point with coordinates (2/3, 0) and no other critical point in the affme plane. 
We note that any function of the form V = f(x, y) n , with n a natural number, 
provides an analytic inverse integrating factor for the system in the whole plane 
M. 2 . Thus, we cannot define the vanishing multiplicity of an inverse integrating 
factor on T. 

We observe that for any natural number n, there exist perturbations of system 
(I22p with at least n limit cycles which bifurcate from the considered homoclinic 
loop T, that is, there is no finite upper bound for the cyclicity of T. In the following 
paragraph, we illustrate this fact by exhibiting a suitable perturbation. We remark 
that the fact that the vanishing multiplicity of an inverse integrating factor on T is 
not defined is coherent with the nonexistence of an upper bound for the cyclicity 

of r. 

For instance, let us fix a natural value n and let e be a nonzero real number 
with |e| small enough. We take e R, i = 1, 2, . . . , n, such that a% ^ % if i ^ j 
and with < a^e < 4/27. For each % we have that the algebraic curve fi = 
with fi := f + diE has an oval with the point (2/3, 0) in its inner region and the 
point (0, 0) in its outer region. Moreover, if e < aj e, the oval defined by fi = 
contains the oval fj = in its inner region. The following system: 

n 

x = —2y, y = — 2x + 3x 2 + ey Y\(f + a i £ ) 

i=l 

is a perturbation of system (1221) and exhibits each one of the ovals defined by the 
curves fi = as periodic orbits. Easy computations show that, if |e| > is small 
enough, the only critical points of the perturbed system are (0,0), which is a hy- 
perbolic saddle point, and (2/3, 0) which is a strong focus. It is also easy to check 
that each oval described by fi = is a hyperbolic limit cycle of the perturbed sys- 
tem which bifurcates from T. We have, in this way, that the considered perturbed 
system has at least n limit cycles which bifurcate from T. 

Acknowledgements. We would like to thank Prof. Cristobal Garcia from 
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